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Transmission  coefficient  measurement  and  improved  sublayer 
material  property  determination  for  thick  underwater  acoustic 
panels:  A  generalization  and  improvement  of  the  ONION  method 

JeanC.  Piquette 

Sava!  Research  Laboratory.  Underwater  Sound  Reference  Detachment,  P.O.  Box  56SJ37,  Orlando,  Florida 
32856-8337 

( Received  1 7  October  1990;  accepted  for  publication  8  April  1992) 

Modifications  of  the  ONION  panel-measurement  method  [  J.  C.  Piquette,  J.  Acoust.  Soc.  Am, 

85,  1029-1040  ( 1989)  ]  that  allow  for  simultaneous  analysis  of  transmitted-  and  reflected-wave 
data  are  described.  The  revised  algorithm  determines  more  reliable  values  for  the  sound  speed 
and  loss  of  the  material  of  each  panel  sublayer  than  does  the  algorithm  that  is  based  exclusively 
upon  analysis  of  the  reflected  wave.  Included  in  the  revised  method  is  a  Taylor  series  expansion 
of  the  sound-speed  function  of  each  layer  about  the  steady-state  driving  angular  frequency. 

This  Taylor  series  is  similar  to  that  used  for  the  loss  function  in  the  original  ONION  method, 
and  is  introduced  here  to  more  accurately  accommodate  the  frequency  variation  of  the  phase 
speed  than  does  the  frequency-independent  model  used  previously.  Descriptions  of  successful 
applications  of  the  revised  ONION  method  to  experimental  data  are  provided.  The  version  of 
the  ONION  method  described  in  this  report  has  recently  been  adopted  as  the  standard  panel 
measurement  method  at  the  Underwater  Sound  Reference  Detachment  of  the  Naval  Research 
Laboratory  (NRL-USRD)  in  Orlando,  FL  for  tests  conducted  in  the  1-  to  20-kHz  frequency 
interval. 

PACS  numbers:  43.20.Fn.  43.20.Px,  43.30.Sf,  43.60.Gk 


INTRODUCTION 

The  ONION  method  is  a  panel  measurement  technique 
that  IS  based  on  least-squares  fitting  of  digitallv  acquired 
transient  pulsed-waveform  data  to  a  multiple-layer  panel 
model.'  '  Readers  who  are  unfamiliar  with  the  technique 
should  consult  the  references  for  complete  descriptions  and 
explanations. 

All  previous  descriptions  of  applications  of  the  method 
have  involved  considerations  of  only  reflected-wave  data. 
The  reason  that  the  transmitted  wave  has  not  heretofore 
been  included  in  the  analysis  is  that  the  problems  involved  in 
accommodating  the  transmitted  wave  are  far  more  formida¬ 
ble  than  are  those  involved  in  accommodating  the  reflected 
wave,  fhe  primary  difficulty  associated  with  the  evaluation 
of  the  transmitted  wave  derives  from  the  very  low  sound 
speeds  that  are  characteristic  of  the  samples  of  interest. 
( These  sound  speeds  can  be  Ic'tS  than  the  sound  speed  in  air. ) 
Thus  (inly  a  very  short  portion  of  the  transmitted  wave  is 
available  for  analysis  prior  to  the  reception  of  unwanted  in¬ 
terfering  waves  (such  as  those  originating  at  the  sample 
edges  or  those  associated  with  measurement  facility  wall 
echo)  In  addition,  since  pulsed  waveforms  are  used  in  the 
tests  of  interest,  the  reflected  waveform  can  be  resolved  into 
convenient  analysis  epochs,  in  which  each  layer  sequentially 
contribiiies  to  the  overall  reflected  waveform.  No  such  cori- 
veiiKii:  suhdiviston  of  the  transmitted  waveform  is  possible, 
since  ven  the  earliest  portions  of  the  transmitted  wave  have 
been  influenced  by  every  layer  of  the  panel. 

Ill  the  present  article,  an  extension  of  the  method  to 


include  analysis  of  transmitted-wave  data  is  described.  The 
extended  ONION-method  algorithm  involves  simultaneous 
least-squares  fitting  of  reflected-wave  data  and  transmitted- 
wave  data  to  a  theoretical  panel  model.  That  is.  the  least- 
squares  calculation  which  is  performed  within  the  extended 
ONION  method  determines  material  properties  for  each  of 
the  panel  sublayers  that  simultaneously  minimize  the  mean- 
squared  error  of  fit  between  the  model  and  the  data  for  both 
the  reflected  waveform  and  the  transmitted  waveform. 

In  the  previous  applications'  '  of  the  ONION  method 
to  only  reflected-wave  data,  the  material  properties  so  deter¬ 
mined  were  considered  to  be  merely  “curve-fitting"  proper¬ 
ties,  and  not  necessarily  the  true  material  properties  of  the 
sublayers.  The  material  properties  determined  in  the  reflect¬ 
ed-wave  analysis  can  often  be  unreliable,  especially  for  layers 
deep  within  the  panel,  due  to  the  fact  that  the  amount  of 
influence  that  a  particular  layer  may  have  on  the  reflected 
wave  can  be  quite  small.  Hence,  any  particular  layer  might 
assume  a  wide  range  of  properties  and  yet  a  good  least- 
squares  minimization  might  nonetheless  be  achieved.  How¬ 
ever,  a  transmitted-wave  analysis  must  necessarily  involve 
information  for  every  sublayer  of  the  panel.  Thus  material 
properties  that  simultaneously  fit  reflected-  and  transmit¬ 
ted-wave  data  are  expected  to  be  more  reliable  than  those 
that  fit  the  reflected  wave  only. 

In  Sec.  I  are  described  certain  improvements  to  the  raw- 
data  handling  algorithms.  In  particular,  this  section  de¬ 
scribes  improvements  to  the  waveform  start-point  determin¬ 
ation  algorithms.  In  Sec.  1!.  modifications  that  have  been 
made  to  the  theoretical  model  incorporated  in  the  ONION- 


468 


J  Acnijst  Soc  Am  92(’)  July  1992 


468 


method  software  are  discussed.  Experiments  conducted  to 
investigate  the  revised  ONION  method  are  discussed  in  Sec. 
III.  A  discussion  of  why  the  material  properties  determined 
by  th'*  revised  method  are  believed  to  be  more  reliable  than 
those  determined  by  the  original  ONION  method  is  given  in 
Sec.  IV.  Section  V  gives  a  summary  and  the  conclusions. 

I.  WAVEFORM  START-POINT  DETERMINATION 

The  theoretical  panel  model  used  to  evaluate  the  reflect¬ 
ed  waveform  in  theONION-method  algorithm  assumes  that 
the  reflected-wave  pulse  has  been  measured  at  the  interface 
between  the  water  medium  and  the  panel  layer  located  clos¬ 
est  to  the  source  of  the  interrogating  wave.  Similarly,  the 
theoretical  panel  model  used  to  evaluate  the  transmitted 
waveform  assumes  that  the  transmitted-wave  pulse  has  been 
measured  at  the  interface  between  the  water  medium  and  the 
panel  layer  farthest  from  the  acoustic  source.  (The  last  panel 
layer  is  usually  a  steel  backing  plate.)  Since  any  practical 
measurement  obviously  requires  an  offset  between  the  rel¬ 
evant  interface  and  the  detecting  hydrophone,  the  experi¬ 
mental  waveforms  so  acquired  must  be  preprocessed,  prior 
to  evaluation  by  the  ONION-method  software.  We  consider 
the  preprocessing  of  each  waveform  separately. 

A.  Reflected  waveform 

As  described  in  Ref  1,  the  reflected  waveform  is  experi¬ 
mentally  determined  using  two  different  measurements.  One 
measurement  involves  capturing  data  in  the  reflection  region 
of  the  panel  with  the  hydrophone  situated  at  a  5-cm  offset 
distance  from  the  relevant  fluid-panel  interface.  The  wave¬ 
form  so  acquired  is  termed  a  "total”  waveform,  since  it  is 
actually  the  sum  of  two  waveforms,  viz.,  the  incident  wave¬ 
form  plus  the  reflected  waveform.  A  second  measurement 
involves  again  capturing  a  waveform  with  the  detecting  hy¬ 
drophone  at  the  same  position  as  was  used  in  acquiring  the 
total  waveform,  but  in  this  second  measurement  the  test  pan¬ 
el  is  removed.  The  waveform  acquired  by  this  measurement 
is  termed  the  "incident"  waveform.  The  reflected  waveform 
is  then  determined  by  digitally  subtracting,  point-by-point, 
the  incident  waveform  from  the  total  waveform. 

The  reflected  waveform  obtained  in  this  way  can  be  con¬ 
sidered  to  consist  of  two  distinct  time  regions.  One  region 
may  be  considered  to  be  a  "quasinull"  region,  resulting  from 
the  point-by-point  subtraction  of  the  incident  waveform 
from  the  total  waveform  in  the  time  region  of  these  wave¬ 
forms  that  is  associated  with  the  round-trip  travel  of  sound 
in  water  between  the  detecting  hydrophone  and  the  relevant 
fluid-panel  interface.  The  second  time  region  represents  the 
actual  reflected  waveform. 

The  incident  and  total  waveforms  in  the  quasinull  time 
region  should  ideally  be  identical,  and  thus,  should  subtract 
to  a  perfect  null.  However,  due  to  the  unavoidable  presence 
of  system  noise,  and  due  to  experimental  difficulties  in  repo¬ 
sitioning  the  detecting  hydrophone  at  the  same  location  in 
both  measurements,  a  residual  nonzero  difference  remains. 
The  preprocessing  software  used  in  the  analysis  of  tin.  earlier 
reports'  ‘  determines  how  much  of  a  time  shift  is  required  to 
make  the  reflected  waveform  appear  to  have  been  measured 


at  the  relevant  interface.  This  is  done  by  (i)  performing  a 
statistical  analysis  of  the  noise  in  the  quasinull  region  and, 
(ii)  comparison  of  this  noisy  quasinull  portion  of  the  wave¬ 
form  to  the  initial  nonzero  portion  of  the  waveform  to  deter¬ 
mine  the  true  start  time  of  the  reflected  wave.  A  time  shift  is 
then  applied  to  the  waveform  to  eliminate  the  quasinull  re¬ 
gion.  Of  course,  this  process  is  necessarily  imperfect,  and 
errors  in  start-point  determination  of  5-10  data  points  (at  a 
4-MHz  data  measurement  rate)  using  this  approach  are  not 
uncommon. 

In  order  to  improve  the  reflected  waveform  start-point 
determination,  a  second  level  of  preprocessing  of  the  data  is 
now  performed.  That  is,  a  second  level  of  preprocessing  that 
takes  as  input  the  already  shifted  waveform  produced  by  the 
first  level  of  preprocessing  described  above,  and  further  im¬ 
proves  the  determination  of  the  start  point  of  the  reflected 
waveform,  is  performed  in  the  revised  version  of  the  soft¬ 
ware.  This  second-level  preprocessing  involves  relying  upon 
the  accuracy  of  the  phase  1  portion  of  the  original  ONION- 
method  algorithm.  (See  Refs.  1  and  2.) 

In  summary,  the  phase  1  portion  of  the  algorithm  is  the 
layer-stripping  (or  "onion-peeling”)  portion.  It  involves  the 
use  of  portions  of  the  reflected  waveform  that  do  not  include 
multiple  internal  reflections,  so  that  the  simple  theoretical 
expressions  for  the  reflection  and  transmission  coefficients 
of  two  semi-infinite  half-spaces  in  intimate  contact  are  appli¬ 
cable.  The  new  waveform  start-point  algorithm  involves,  (i) 
attempting  a  candidate  time  shift  of  the  experimental  reflect¬ 
ed  waveform  corresponding  to  a  predetermined  discrete 
number  of  data  points,  (ii)  applying  the  phase- 1  layer-strip- 
ping  algorithm  to  the  shifted  reflected  waveform  to  deduce 
an  approximate  sound  speed  and  loss  for  the  panel  layer 
situated  closest  to  the  acoustic  source  and.  ( iii )  computing  a 
theoretical  pulsed  waveform,  based  on  the  stripped  values  of 
the  material  properties.  The  root-mean-squared  error  of  fit 
between  the  model  waveform  and  the  shifted  experimental 
waveform,  evaluated  up  to  that  point  in  time  at  which  the 
first  internal  reflection  is  expected  to  occur,  is  then  comput¬ 
ed  and  stored.  Next,  another  candidate  time  shift  is  applied 
to  the  experimental  reflected  waveform,  and  the  above-de¬ 
scribed  calculations  are  run  for  the  new  shift.  The  second- 
level  preprocessing  algorithm  proceeds  in  this  manner  until 
it  has  evaluated  and  stored  root-mean-squared  errors-of-fit 
that  correspond  to  a  sequence  of  candidate  shifts  that  vary  in 
discrete  amounts  from  1/4  cycle  (of  the  steady-state  driving 
frequency)  earlier  in  time,  up  to  1/4  cycle  later  in  time,  rela¬ 
tive  to  the  start  point  of  the  waveform  produced  by  the  first- 
level  quasinull  elimination  calculation  described  above.  The 
shift  that  produces  the  least  root-mean-squared  error  of  fit  is 
taken  to  be  the  correct  shift,  and  this  shift  is  then  applied  to 
the  experimental  reflected  waveform  and  retained.  The  en¬ 
tire  waveform-shifting  calculation  described  here  requires 
approximately  1-CPU  min  to  perform  on  a  DEC  Micro- 
Vax™  3XXX  series  workstation  computer. 

B.  Transmitted  waveform 

The  waveform-shifting  technique  described  above  for 
the  reflected  waveform  is  not  applicable  to  the  transmitted 
waveform.  This  is  due  to  the  fact,  previously  mentioned,  that 
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the  onion-peeling  portion  of  the  algorithm  does  not  produce 
reliable  properties  for  layers  situated  deep  within  the  panel. 
(Unlike  the  reflected  waveform,  whose  initial  data  points 
depend  only  upon  properties  of  the  first  panel  layer,  all  of  the 
data  points  of  the  transmitted  waveform  depend  upon  the 
properties  of  all  panel  layers.)  Thus  the  time  shift  of  the 
transmitted  waveform  must  be  effected  by  a  different,  and 
less  accurate,  method. 

An  additional  difficulty  associated  with  the  determina¬ 
tion  of  the  time  location  of  the  start  point  of  the  transmitted 
wave  is  caused  by  the  ( unknown )  time  delay  associated  with 
the  time  of  flight  of  the  transmitted  '^ave  from  the  front 
fluid-panel  interface  (i.e.,  the  interface  closest  to  the  acous¬ 
tic  source)  to  the  back  fluid-panel  interface  (i.e.,  the  inter¬ 
face  farthest  from  the  acoustic  source).  This  delay  is  experi¬ 
mentally  determined  by  performing  two  waveform 
measurements  in  the  transmission  region  of  the  panel.  We 
discuss  the  measurement  process  with  the  aid  of  Fig.  1. 

In  Fig.  1  are  presented  the  two  waveforms  which  must 
be  measured  in  order  to  make  the  transmitted  wave  start- 
point  determination  presently  under  consideration.  The  first 
displayed  waveform  represents  the  incident  wave  as  it  ap¬ 
pears  in  the  transmission  region  with,  the  sample  panel  re- 
moved.  Also  shown  is  the  transmitted  waveform  as  it  is  mea¬ 
sured  with  the  sample  panel  in  place.  (Note  that  the 
waveforms  are  sketched  to  suggest  the  dispersive  nature  of 
the  sample  material.  In  particular,  note  the  different  slopes 
depicted  for  the  initial  nonzero  values  of  the  incident  wave¬ 
form  as  compared  to  those  of  the  transmitted  waveform.) 
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interrogating  wave).  Ar  -  digili/ed  sample  spacing  in  the  time  domain 
( tvpicallv  t)  2?  /fs ).  n  ,,  , ,  ,,,  data  point  number  fif  the  start  point  of  the 
nil. idem  waveform  as  measured  in  the  transmission  region  with  the  viniple 
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The  waveforms  are  displayed  with  coincident  time  coordi¬ 
nates,  with  t  =  0  defined  to  be  the  arrival  time  of  the  incident 
wave  at  the  “front”  panel  interface;  i.e.,  the  sample  surface 
located  closest  to  the  source  of  the  interrogating  wave.  Each 
waveform  has  an  associated  null  region,  which  is  subdivided 
in  the  figure  into  two  subregions.  The  time  region  is  the 
time  of  flight  of  the  incident  wave  in  passing  through  the 
region  of  water  that  is  located  where  the  panel  was  prior  to 
its  removal.  If  we  let  /  denote  the  overall  sample  thickness, 
and  let  c„.  denote  the  speed  of  sound  in  water,  then 
■  (The  value  of  /  is  determined  in  the  ONION 
method  with  the  help  of  an  underwater  video  camera. ) 

The  /;,  time  region  denotes  the  time  of  flight  of  the  inci¬ 
dent  wave  from  the  “back”  panel  interface  ( i.e.,  the  interface 
farthest  from  the  acoustic  source)  to  the  detecting  hydro- 
piione  in  the  transmission  region.  Due  to  the  uncertainty  in 
the  location  of  the  effective  acoustic  center  of  the  hydro¬ 
phone,  this  time  cannot  be  reliably  determined  using  a  dis¬ 
tance  measurement  of  the  sort  discussed  above  to  determine 
the  quantity  However,  the  quantity  can  be  determined 
indirectly,  as  we  will  see  presently. 

For  the  transmitted  waveform,  the  quantity  ti,.,,,..,  de¬ 
notes  the  ( unknown )  time  of  flight  of  the  transmitted  wave 
from  the  front  panel  interface  to  the  back  pane!  interface, 
through  the  actual  material  of  the  sample.  The  quantity 
has  the  same  meaning  as  discussed  above  in  connection  with 
the  incident  waveform. 

Note  from  Fig.  1  that,  in  order  to  make  the  transmitted 
waveform  appear  to  have  been  measured  at  the  back  panel 
interface,  the  waveform  must  be  shifted  (to  the  “left”)  by  a 
data-point  number  that  is  equivalent  to  the  time  l,,.  This 
quantity  is  determined  from  the  incident-wave  measurement 
in  the  following  way:  First,  the  start  point  of  the  incident 
waveform,  as  measured  in  the  transmission  region,  is  deter¬ 
mined  by  performing  a  similar  type  of  statistical  comparison 
of  the  null  part  of  the  waveform  with  the  initial  nonzero 
portion  of  the  waveform  as  is  used  to  make  the  first-level 
start-point  determination  for  the  reflected  wave.  (Of  course, 
this  method  is  significantly  more  reliable  when  applied  to  the 
incident  wave  than  when  it  is  applied  to  either  the  reflected 
or  transmitted  wave  because  of  the  much  ““eater  signal-to- 
noise  ratio  of  the  incident  wave.)  The  d.Ua-point  number 
that  specifies  the  start  point  of  the  iu'ident  waveform  as 
measured  in  the  transmission  regio'^  ith  the  panel  removed 
is  termed  in  Fig.  1 .  The  quantity  At  of  Fig.  1  denotes 

the  sample  spacing  in  the  time  domain  of  the  digitized  wave¬ 
forms.  (Typically,  At  =  0.25  ps  in  the  ONION  method.) 
The  quantity  t,,  can  be  determined  using  the  symbols  defined 
here  by  substitution  int»  the  formula  t,,  =  At  —  t,,^ 

( see  Fig.  1 ) ,  where,  discussed  above,  t„ ,,  =  I  /c„ .  Thus  the 
required  data-poi  ii  shift"* '  of  the  transmitted  wave  is  given 
in  terms  of  known  quantities  by  the  expression  t,,/At,  which 
is  equivalent  lo  -  ^,,,/At,  or  -  / /c„  At. 

II.  THE  CAUSALITY  PROBLEM  AND  MODEL 
IMPROVEMENTS 

As  was  discussed  in  a  previous  article.’  the  ONION- 
method  model  is  noncausal.  This  difficulty  arises  from  the 
treatment  of  the  frequency  dependence  of  the  sound  speed 
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and  loss  functions  in  the  model.  The  model  treats  the  loss 
function  of  the  panel  materials  as  frequency  dependent,  but 
treats  the  sound  speed  as  frequency  independent.  This  treat¬ 
ment  of  these  functions  is  inconsistent  with  the  requirements 
of  causality.'’  Although  provisions  to  account  for  the  non¬ 
causality  of  the  model  are  present  in  the  software  that  imple¬ 
ments  the  ONION-method  algorithm,  it  is  nonetheless 
worthwhile  to  investigate  modifications  of  the  model  that 
can  render  it  consistent  with  the  causality  principle.  This  is 
particularly  important  in  the  present  work  since  it  is  almost 
certainly  not  sufficient  to  merely  deduce  “curve-fitting” 
properties  if  the  simultaneous  fitting  of  reflection  and  trans¬ 
mission  data  is  desired,  as  is  the  case  here. 

In  order  to  investigate  how  the  GNION  model  might  be 
rendered  consistent  with  the  cau.sality  principle,  a  number  of 
attenuation-dispersion  pairs  available  in  the  geophysical  li¬ 
terature"  ''  were  considered.  The  software  that  implements 
the  model  calculations  was  modified  to  successively  incorpo¬ 
rate  each  of  these  literature  models.  The  modified  software 
was  then  used  to  analyze  panel  calibration  data.  Unfortu¬ 
nately,  none  of  these  models  proved  to  be  any  more  success¬ 
ful  in  fitting  the  data  than  the  noncausal  model  originally 
.o.isidered.'  That  is,  the  mean-squared  error  of  fit  between 
the  model  and  the  data  was  not  found  to  be  improved  by 
using  these  models  from  the  literature. 

This  result  is  not  as  surprising  as  it  might,  at  first,  ap¬ 
pear.  Although  the  literature  models  satisfy  the  causality 
principle  e.xactly,  this  e.xact  behavior  is  achieved  at  the  cost 
ofassummg  a  specific  (and  quite  simple)  functional  form  for 
the  loss  function.  (For  e.xample,  the  model  considered  by 
Futterman' '  assumes  the  loss  function  varies  linearly  with 
frequency. )  The  sound-speed  variation  with  frequency  is 
then  deduced  using  the  appropriate  Hilbert  transform."  This 
Hilbert  transform  calculation  assures  that  the  causality  prin¬ 
ciple  will  be  satisfied  exactly.  However,  if  the  experimental 
data  exhibit  behavior  that  is  inconsistent  with  the  assumed 
frequency  dependence  of  the  loss  function,  the  model  will 
nonetheless  fail  to  fit  the  data  well. 

Of  course,  considering  the  manner  in  which  typical  pan¬ 
el  sublayers  are  fabricated,  subsuming  such  complexities  as 
voids  and  seams,  there  is  no  compelling  reason  to  suppose 
that  a  linear  variation,  or  any  other  simple  functional  vari¬ 
ation,  of  the  loss  function  with  frequency  will  be  satisfac¬ 
tory.'"  In  fact,  considering  the  unknowns  involved  in  sub¬ 
layer  fabrication,  such  as  void  size,  shape,  and  distribution,  it 
IS  unlikely  that  a  simple  loss  function  exists  that  can  accu¬ 
rately  accommodate  all  cases  of  current,  or  potential  future, 
interest. 

The  ONION-method  model,  despite  being  noncausal, 
can  accommodate  a  w  ide  variety  of  materials  due  to  the  fact 
that  no  strict  assumptions  are  made  in  the  functional  form  of 
the  loss  function.  That  is.  the  loss  function  is  modeled  by  a 
Fay  lor  series,  with  unknown  expansion  coefficients,  truncat¬ 
ed  at  the  quadratic  term.  This  model  is,  of  course,  compati¬ 
ble  with  any  data  that  exhibit  constant,  linear,  or  quadratic 
dependence  of  the  loss  function  on  frequency.  However,  the 
model  IS  also  at  least  approximately  consistent  with  much 
more  complicated  frequency  variations,  provided  that  the 
quadratic-term  truncation  of  the  Tay  lor  series  yields  a  rea¬ 


sonable  approximation  to  the  true  loss  function.  [Since  the 
frequency  content  of  the  interrogating  pulses  used  in  the  ex¬ 
periments  performed  to  implemciu  ONION-method  mea¬ 
surements  covers  a  relatively  narrow  0-  to  40-kHz  band 
(with  no  significant  information  below  about  1(X)  Hz),  and 
since  the  Taylor  series  is  expanded  about  the  steady-state 
interrogating-wave  frequency,  this  truncated  Taylor  series  is 
likely  to  be  a  reasonable  approximation  to  a  wide  class  of  loss 
functions. )  The  ability  of  the  model  to  accommodate  a  wide 
variety  of  loss  function  frequency  variations  serves  to  ex¬ 
plain  its  success  in  evaluating  panel  calibration  data,  despite 
the  inconsistency  with  the  causality  principle. 

In  view  of  the  superior  performance  of  the  ONION- 
method  model  compared  to  the  literature  models  in  evaluat¬ 
ing  panel-test  data,  and  considering  the  fact  that  no  compell¬ 
ing  dispersion-attenuation  pair  is  available  for  the  complex 
materials  used  in  panels,  it  was  decided  to  abandon  at  the 
present  time  the  attempt  to  render  the  ONION  model  strict¬ 
ly  consistent  with  the  causality  principle.  It  was  decided  in¬ 
stead  to  introduce  additional  flexibility  into  the  model  to  at 
least  permit  causal  behavior.  This  is  accomplished  in  the 
present  version  of  the  ONION-method  model  by  introduc¬ 
ing  a  truncated  Taylor  series  for  the  frequency  variation  of 
the  sound-speed  function  that  is  similar  to  that  introduced 
previously  for  the  loss  function. 

The  Taylor  series  previously  introduced  to  accommo¬ 
date  the  frequency  variation  of  the  loss  function  has  the 
form' 
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and  the  additional  Taylor  series  hereby  introduced  to  better 
accommodate  the  frequency  variation  of  the  phase  speed 
function  is 
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In  these  equations,  w  is  the  (variable)  angular  frequency,  w„ 
is  the  (constant)  steady-state  driving  angular  frequency  of 
the  interrogating  wave,  A  is  the  wavelength  of  the  acous¬ 
tic  wave  in  the  layer  of  interest  at  angular  frequency  w,,,  c„  is 
the  sound  speed  in  water,  and  (Wo  )  is  the  .sound  speed  in 
the  layer  of  interest  at  w,,.  Theparameter  w  isan  integer  used 
to  number  each  layer.  The  Taylor  series  expansion  constants 
air'- a'"'-  ).  ^  'i'"’.  and  b  are  the  parameters 

that  are  adjusted  using  a  least-squares  fitting  process  that  is 
similar  to  that  described  in  Refs.  1  and  2. 

In  Fig.  2  is  presented  a  block  diagram  that  summarizes 
the  revised  ONION  method.  The  method  requires  two 
phases.  As  with  the  earlier  implementation,'  ’  phase  1  of  the 
present  implementation  uses  portions  of  the  reflected  time 
waveform  that  are  free  from  multiple  internal  reflections,  so 
that  the  simple  theoretical  expressions  for  the  reflection  and 
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[  Ki,  ’  Block  diajiruin  ol  revised  ONION-mcihod  algorithm.  As  with  the 
earlier  implementation  i  see  Rel's.  1  aiul  ’ ),  the  process  requires  iwopliases. 
Phase  1  represents  ihc  "laver-stripping''  ( or  "onion-peeling"  )  phase,  dur¬ 
ing  which  data  that  are  free  from  multqile  mternal  reflections  are  analyzed, 
using  the  simple  e.’ipressions  for  the  reflection  and  transmission  coefficients 
for  two  senii-intinite  media.  This  analysis  produces  appro.ximaie  sound 
speeds  and  losses  for  each  l.iycr.  In  phase  2.  a  nonlinear  least-squares  calcu¬ 
lation  IS  used  Iteratively  to  improve  the  appro.vimate  paranieter  values.  Dur¬ 
ing  this  phase,  a  simultatieous  fit  of  experimental  reflected  and  transmitted 
waveforms  to  appropriate  model  waveforms  is  performed  I'nlikelhe  meth¬ 
od  of  Refs  I  and  2,  the  model  used  in  the  present  implementation  includes  a 
l  ay  lor  series  e.xpaiision  for  the  phase-speed  function  for  each  layer  in  addi- 
t.on  to  that  previously  used  for  the  loss  function.  Here.  Ain  is  the e.xperi- 
niental  piilsed-incident  time  waveform.  A,  I  r)  is  the  experimental  puKed- 
reflected  tune  w.iveform.  and  A,  1  /)  is  the  experimental  pulsed-transmitted 
lime  waveform  Similarly,  P^tl)  is  the  computed  pulsed-reflected  lime 
w.iveform  aiul  A,  i  n  is  the  computed  pulsed-lransmilled  lime  waveform, 
based  on  applying  the  transfer  funciion  of  mulliple-layer  theory,  suitably 
modified  to  incorporate  the  two  ahove-meiilioned  Tay  lor  senes,  to  Air) 


transmissKin  coefficients  for  two  semi-infinite  media  in  inti¬ 
mate  cyintact  are  applicable.  This  calculation  is  used  to  de¬ 
duce  appro.ximate  starting  parameter  values.  In  phase  2  a 
nonlinear  least-squares  fitting  algorithm  is  used  iteratively 
to  produce  improved  parameter  values.  This  process  pro¬ 
duces  parameter  values  that  are  most  consistent  (in  a  least- 
squares  sense )  w  ith  the  available  data.  In  view  of  the  dispa¬ 
rate  nature  of  the  two  sets  of  waveforms  being  fitted  (i.e.,  the 
reflected  and  transmitted  waveforms),  the  sound  speeds  and 
losses  obtained  in  this  way  are  e.xpected  to  be  more  reliable 
estimatec  of  the  true  effective  sound  speeds  and  losses  of  each 
layer  than  are  those  deduced  previously'  by  only  fitting  rc- 
flected-wave  data.  This  will  be  discussed  further  in  See,  IV. 

It  IS  important  at  this  point  to  mention  that  in  imple¬ 
menting  the  simultaneous  least-squares  minimization  pro¬ 
cess  of  the  phase-2  portion  of  the  algorithm,  a  weighting 


scheme  is  used  to  put  the  two  fitted  waveforms  on  an  equal 
footing.  That  is.  since  the  amplitude  of  the  transmitted  wave 
is  often  much  lower  than  that  of  the  reflected  wave  (some¬ 
times  more  than  40  dB  lower),  if  no  account  were  taken  of 
this  lower  amplitude  then  the  differences  between  the  model 
and  the  data  for  the  transmitted  wave  would  be  completely 
ineffectual  in  driving  the  calculations  of  model  corrections. 
To  correct  for  this,  the  ratio  of  the  mean-squared  amplitude 
of  the  reflected  wave  to  that  of  the  transmitted  wave  is  used 
as  a  weighting  factor  to  put  the  two  waves  on  an  equal  foot¬ 
ing.  This  weighting  factor  is  used  in  evaluating  the  contribu¬ 
tions  of  the  transmitted-wave  data  in  the  least-squares  fitting 
process. ' ' 


III.  EXPERIMENTS 

Experimental  measurements  were  undertaken  to  inves¬ 
tigate  the  effectiveness  of  the  revised  ONION  method.  As 
was  done  in  the  earlier  reports.'  '  tests  were  carried  out  on 
two  different  samples.  One  sample  contains  simple  homoge¬ 
neous  layers,  while  the  other  contains  complex  layers  fabri¬ 
cated  from  voided  rubbers  affixed  to  a  steel  backing  plate. 
Measurements  performed  on  the  sample  containing  simple 
homogeneous  layers  covered  the  frequency  range  3-10  kHz. 
Measurements  performed  on  the  sample  containing  voided- 
rubber  sublayers  covered  the  frequency  range  1-10  kHz.  In 
both  cases,  the  source-to-sample  separation  used  was  200 
cm;  sample-to-detector  separations  were  5  cm  from  each  rel¬ 
evant  interface  on  both  the  reflected  and  transmitted  sides  of 
the  sample.  Transmitted  waveforms  were  acquired  using  a 
configuration  obtained  by  rotating  the  sample  180°  with  re¬ 
spect  to  the  configuration  used  to  obtain  the  reflected  wave¬ 
forms.  As  has  been  observed  previously."*  this  configuration 
(in  which  the  backing  plate  faces  the  acoustic  source)  is  a 
better  experimental  realization  of  infinite-sample  transmit- 
led-wave  theory  than  is  the  standard  configuration  (in 
which  the  backing  plate  is  located  on  the  side  of  the  panel 
opposite  the  acoustic  source).  That  is.  fitting  errors  are 
found  to  be  less  using  transmitted  waveforms  acquired  with 
the  backing  plate  facing  the  source  than  when  using  trans¬ 
mitted  waveforms  acquired  when  the  backing  plate  is  oppo¬ 
site  the  source. 

The  panel  containing  simple  homogeneous  layers  con¬ 
sists  of  one  layer  of  polymethylmethacrylate  of  2.54-cm  ( l- 
in.)  thickness,  followed  by  a  water  layer  of  2.54-cm  thick¬ 
ness  and  a  steel  layer  of  0.95-cm  (3/8-in.)  thickness.  The 
sample  containing  voided-rubber  layers  consists  of  a  layer  of 
density  0.78-g/cm '  and  4.84-cm  thickness  which  is  laminat¬ 
ed  onto  a  second  layer  of  density  0.91-g/cm'  and  4.84-cm 
thickness,  followed  by  a  third  layer  of  density  0.&2-g/cm' 
and  2.84-cm  thickness.  The  third  layer  is  laminated  onto  a 
standard  steel  support  plate  of  ().95-cm  (3/8-in.)  thickness. 
BtUh  samples  have  square  lateral  area.  76  cm  (30  in.)  on  a 
side. 

In  F'igs.  3  and  4  are  presented  representative  waveforms 
for  each  of  these  samples.  In  Fig.  3  are  presented  representa¬ 
tive  waveforms  for  the  sample  containing  simple  homoge¬ 
neous  layers.  Test  frequency  ts  8  kHz,  In  Fig,  4  are  presented 
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F’lG.  .V  Typical  cxpcrmicMUai  and  model  waveforms  for  the  sample  containing  three  simple  homogeneous  layers.  Frequency  is  8  kHz.  Solid  line — experimen¬ 
tal  waveforms.  Dashed  line-  numencally  computed  waveforms  based  on  least-squares  fitting  to  the  experimental  data  available  in  the  indicated  data  win¬ 
dows.  Vertical  lines  delineate  data  windows,  (a)  Transmitted  and  tb)  reflected. 
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representative  waveforms  for  the  sample  containing  voided- 
rubber  sublayers.  Test  frequency  is  2  kHz. 

As  can  be  seen,  the  fit  of  the  model  waveforms  to  the 
data  within  the  indicated  data  windows  is  excellent  in  both 
cases.  The  model  waveforms  and  experimental  waveforms 
for  the  sample  containing  simple  homogeneous  layers  can  be 
compared  outside  the  data  windows  to  deduce  that  edge  dif¬ 
fraction  is  of  minor  importance  in  the  transmitted  wave 
[note  Fig.  3(a)]  but  is  somewhat  more  significant  in  the 
reflected  wave  [note  Fig.  3(b)].  [The  same  type  of  compari¬ 
son  cannot  be  done  for  the  waveforms  acquired  from  the 
samples  containing  voided  viscoelastic  sublayers  (i.e.,  the 
waveforms  of  Fig.  4),  due  to  the  need  to  use  all  of  the  avail¬ 
able  data  in  the  fitting  process.  This  is  necessary  due  to  the 
very  low  sound  speeds,  and  greater  layer  thicknesses,  of  this 
sample.  ]  Calculations  show  that  the  average  root-mean- 
squared  error  of  fit  for  the  waveforms  of  Fig.  3  within  the 
data  windows  is  5.029f ,  while  that  for  the  waveforms  of  Fig. 
4  is  4.31'%-.  These  fitting  errors  are  typical  of  those  obtained 
for  the  entirety  of  the  data. 

In  Figs.  5  and  6  are  presented  graphs  that  summarize  the 
results  obtained  for  all  of  the  measurements.  That  is,  these 
graphs  present  the  transmission  and  reflection  coefficients  as 
a  function  of  test  frequency  for  each  case.  In  Fig.  5  are  pre¬ 
sented  results  for  the  sample  containing  simple  homoge¬ 


neous  layers.  In  Fig.  6  are  presented  results  for  the  sample 
containing  voided-rubber  sublayers. 

In  Tables  I  and  II  are  presented  the  material  properties 
deduced  for  each  of  these  samples  at  the  frequencies  used  in 
obtaining  the  data  presented  in  Figs.  3  and  4,  respectively.  In 
Table  I  are  presented  results  for  the  sample  containing  sim¬ 
ple  homogeneous  layers.  (The  phase  speed  for  PMM  that  is 
available  in  the  literature,'*'  presented  here  for  comparison, 
is  2.68 X  10'  cm/s.)  In  Table  II  are  presented  results  for  the 
sample  containing  voided-rubber  sublayers.  Note  in  Table  II 
that  the  sound  speeds  decline  with  increasing  layer  number, 
while  losses  increase  with  increasing  layer  number.  Such  be¬ 
havior  is  a  typical  design  target  of  panels  of  this  type. 

IV.  DISCUSSION 

We  next  discuss  the  reasons  it  is  believed  that  the  mate¬ 
rial  properties  deduced  by  the  present  method  are  in  fact 
reasonable  estimates  of  the  true  material  properties  of  the 
panel  layers.  This  point  of  view,  of  course,  is  considerably 
different  from  that  taken  in  the  previous  reports'  '  of  work 
on  the  ONION  method. 

First,  it  is  worthwhile  to  recall  the  reason  that  in  the 
earlier  reports  the  material  properties  obtained  were  only 
considered  to  be  “curve-fitting”  properties,  and  not  neces- 
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earliest  non/ero  values  of  the  transmitted  waveform  have 
been  influenced  by  all  the  sublayers  of  the  panel.  Note 
further  that,  if  a  sonntl  speed  for  one  of  the  layers  were  to 
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wander  very  far  from  its  true  value,  an  unacceptable  phase 
shift  would  be  introduced  into  the  computed  transmitted 
waveform.  That  is,  the  start  time  of  the  computed  transmit¬ 
ted  wave  would  be  incompatible  with  that  of  the  observed 
transmitted  w'ave.  Hence,  a  good  fit  between  model  and  data 
could  not  be  achieved.  Of  course,  a  too-great  phase  speed  for 
one  layer  could  conceivably  be  compensated  by  a  too-low 
phase  speed  for  another  layer.  However,  in  this  case,  it  seems 
unlikely  that  the  resulting  waveshape  of  the  computed  trans¬ 
mitted  wave  could  properly  accommodate  that  of  the  mea¬ 
sured  waveform. 

Confidence  in  the  results  has  also  been  achieved  by  at¬ 
tempting  certain  numerical  tests.  In  one  such  test,  the  layer- 
striping  (or  onion-peeling)  portion  of  the  algorithm  was 
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over-ridden,  and  relatively  arbitrary  starting  values  of  the 
material  properties  were  used.  For  example,  in  one  test,  in¬ 
volving  the  data  acquired  from  the  sample  containing  simple 
homogeneous  layers,  a  phase  speed  for  PMM  that  is  20% 
greater  than  its  true  value,  and  a  phase  speed  for  water  equal 
to  one-half  its  true  value,  were  used  as  initial  model  param¬ 
eters.  Thus  the  fitting  algorithm  had  every  opportunity  to 
wander  into  wildly  inappropriate  regions  of  parameter  space 
in  this  test.  However,  the  algorithm  nonetheless  converged 
to  phase  speeds  for  the  layers  that  were  within  15%  of  the 
true  values.  Losses  were  also  accurately  determined  to  have 
a  negligible  value. 

In  a  second  test  using  these  data,  the  water  layer  was 
analytically  subdivided  into  sublayers,  with  each  such  sub¬ 
layer  having  one-half  the  true  thickness  of  the  layer  in  ques¬ 
tion.  Each  of  these  sublayers  was  given  an  initial  sound  speed 
equal  to  one-half  the  true  sound  speed  of  water,  and  the 
properties  of  each  sublayer  were  permitted  to  be  indepen¬ 
dently  adjusted  by  the  software.  The  final  phase  speeds  de¬ 
duced  for  the  two  water  sublayers  were  within  5%  of  each 
other,  and  were  also  each  within  1 5%  of  the  true  phase  speed 
of  water  (Again,  loss  was  accurately  determined  to  have  a 
negligible  value.)  Thus  the  algorithm  is  seen  to  be  robust 
against  poor  initial  values  for  the  layer  pioperties. 

One  final  test  of  the  accuracy  of  the  properties  deduced 
by  the  revised  ONION  method  will  be  described.  In  recent 
tests  conducted  at  our  laboratory,  a  sample  panel  designed 
for  a  decoupling  application  in  an  acoustic  array  was  evalu¬ 
ated.  The  intended  application  of  the  panel  involves  a  situa¬ 
tion  in  which  a  vibrating  metal  surface  in  contact  with  water, 
but  backed  by  air,  is  desired  to  be  acoustically  isolated.  The 
desired  isolation  is  specified  by  a  velocity-reduction  design 
requirement;  i.e..  the  velocity  amplitude  on  the  decoupler 
surface  must  be  a  specified  amount  smaller  than  that  on  the 
metal  surface.  The  effectiveness  of  the  candidate  material  in 
providing  the  desired  isolation  was  experimentally  deter¬ 
mined  in  two  ways;  ( i )  using  a  special  test  rig  consisting  of  a 
metallic  plate,  backed  by  air,  onto  which  the  sample  is  al- 
fi.xed,  and  ( ii )  performing  a  panel  test  using  the  revised  ON¬ 
ION  method.  In  (i).  the  sample  mat'.-'ial  is  immersed  in 
water,  and  a  direct  observation  of  the  velocity  reduction  is 
made.  In  (li ),  a  plastic  backing  plate'"  is  affixed  to  the  sam¬ 
ple.  and  the  hacking  plate-sample  combination  is  immersed 
in  water.  The  velocity  reduction  for  an  air-hacked  metal 
plate  is  then  deduced  by  analytically  removing  the  plastic 
backing  plate  ami  water  backing  of  test  (ii),  by  using  the 
material  properties  deduced  m  (ii)  by  the  revised  ONION 
software,  and  then  analytically  inserting  an  air-backed  metal 
plate.  The  velocity  reduction  so  obtained  is  found  to  be  in 
reasonably  good  agreement  with  that  obtained  from  the  di¬ 
rect  observation  in  ( i )  .  In  particular,  the  two  velocity-i  educ¬ 
tion  curves  were  found  to  agree  to  within  about  1  dB  over  a 
frequency  interval  exceeding  an  octave.  This  agreement  is 
quite  good  in  view  of  the  fact  that,  due  to  manufacturing 
difficulties,  the  two  tested  samples  were  similar  fat  not  iden¬ 
tical.  Such  agreement  would  be  unlikely  to  occur  if  the  mate¬ 
rial  properties  used  in  these  calculations  were  poor  estimates 
of  the  true  material  properties. 

It  will  he  noted  that  in  the  results  presented  in  Figs.  5 


and  6  no  error  bars  are  given,  unlike  similar  graphs  present¬ 
ed  previously.’  The  reason  for  this  is  that  the  fitting  errors 
deduced  using  the  standard  propagation-of-errer  methods 
previously  used  might  not  actually  provide  reliable  error  es¬ 
timates.  This  is  due  to  the  fact  that  there  are  known  sources 
of  systematic  error  present  in  the  measurements.  This  sys¬ 
tematic  error  is  associated  with  the  rather  significant  edge- 
effects  that  are  known  to  be  present  in  these  measurements, 
especially  in  the  transmitted-wave  case.  ( See  Ref.  18).  These 
edge  effects  are  associated  with  edge  diffraction  and  with 
surface  waves  induced  in  the  sample  surfaces  due  to  the  pres¬ 
ence  of  the  sample  edges.  ( Recall  that  it  is  necessary  to  in¬ 
clude  all  the  available  measured  waveform  in  the  transmit¬ 
ted  wave  analysis  due  to  the  low  sound-speed  and  large  layer 
thicknesses  of  the  voided-rubber  sample.)  The  matter  of 
how  to  properly  account  for  this  systematic  source  of  error  is 
the  subject  of  further  research. 

V.  SUMMARY  AND  CONCLUSIONS 

A  revised  version  of  the  ONION-method  software  that 
simultaneously  fits  reflected-  and  transmitted-wave  data  has 
been  described.  The  revised  version  incorporates  Taylor  se¬ 
ries  expansions  of  both  the  loss  function  and  phase-speed 
function  of  each  panel  layer.  The  modified  ONION-method 
model,  despite  not  being  strictly  causal,  has  been  found  to  fit 
experimental  data  more  accurately  than  several  exactly 
causal  models  available  in  the  literature.  The  material  prop¬ 
erties  deduced  for  each  panel  layer  by  the  revised  method  are 
believed  to  be  reasonably  accurate  determinations  of  these 
nroperiics.  and  are  no  longer  regarded  as  merely  being 
“curve-fitting”  properties,  as  was  the  case  when  only  the 
reflected  waveform  was  used  as  the  basis  of  the  least-squares 
analysis.  It  is  concluded  that  the  revised  ONION  method 
provides  reasonably  accurate  determinations  of  reflection 
and  transmission  coelTicients,  as  well  as  sound  speeds  and 
losses  for  each  panel  layer,  as  a  function  of  temperature, 
pressure,  and  frequency. 

In  closing,  it  is  worthwhile  to  point  out  that  ihe  revised 
version  of  the  ONION  method  described  in  this  report  has 
recently  been  adopted  as  the  standard  panel  measurement 
method  at  the  Underwater  Sound  Reference  Detachment  of 
the  Naval  Research  Laboratory  (NRL-USRD)  in  Orlando, 
FL  for  panel  tests  conducted  in  the  1-  to  20-kHz  frequency 
interval. 
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